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Abstract: Grassland degradation is influenced by climate change and human activities, and has become a 
major obstacle for the development of arid and semi-arid areas, posing a series of environmental and 
socio-economic problems. An in-depth understanding of the inner relations among grassland vegetation 
dynamics, climate change, and human activities is therefore greatly significant for understanding the 
variation in regional environmental conditions and predicting future developmental trends. Based on 
MODIS (moderate resolution imaging spectroradiometer) NDVI (normalized difference vegetation index) 
data from 2000 to 2020, our objective is to investigate the spatiotemporal changes of NDVI in the Xilin 
Gol grassland, Inner Mongolia Autonomous Region, China. Combined with 12 natural factors and human 
activity factors in the same period, the dominant driving factors and their interactions were identified by 
using the geographic detector model, and multiple scenarios were also simulated to forecast the possible 
paths of future NDVI changes in this area. The results showed that: (1) in the past 21 a, vegetation cover 
in the Xilin Gol grassland exhibited an overall increasing trend, and the vegetation restoration (84.53%) 
area surpassed vegetation degradation area (7.43%); (2) precipitation, wind velocity, and livestock number 
were the dominant factors affecting NDVI (the explanatory power of these factors exceeded 0.4). The 
interaction between average annual wind velocity and average annual precipitation, and between average 
annual precipitation and livestock number greatly affected NDVI changes (the explanatory power of these 
factors exceeded 0.7). Moreover, the impact of climate change on NDVI was more significant than human 
activities; and (3) scenario analysis indicated that NDVI in the Xinlin Gol grassland increased under the 
scenarios of reduced wind velocity, increased precipitation, and ecological protection. In contrast, 
vegetation coverage restoration in this area was significantly reduced under the scenarios of unfavorable 
climate conditions and excessive human activities. This study provides a scientific basis for future 
vegetation restoration and management, ecological environmental construction, and sustainable natural 
resoutce utilization in this area. 
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1 Introduction 


Grassland ecosystems cover 30%—40% of the Earth's total land area, contain about 20% of the 
global soil carbon pool (Subramanian et al., 2020), and have a large carbon sink potential (Xin et 
al., 2020). Thus, grassland ecosystems greatly influence the global carbon cycle and regional 
economy (Wang et al., 2019; Yu et al., 2019). In recent decades, grassland degradation has been 
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influenced by climate change and human activities, and has become a major factor restricting 
development in arid and semi-arid areas. Grassland degradation poses a series of environmental 
and socio-economic problems (He et al., 2015; Gu et al., 2018). Dynamic grassland vegetation 
monitoring and corresponding driver analysis have always been an important topic in 
globalization research (Pan et al., 2018; Du et al., 2020). Assessing and analyzing the relative 
roles of climatic factors and human activities in the changes of grassland ecosystem are critical, 
which has important scientific and practical significance for the restoration and protection of 
grassland ecosystems, as well as the construction of regional ecological security barriers. 

Normalized difference vegetation index (NDVI) is the best indicator for monitoring terrestrial 
vegetation status (Alfredo, 2016; Burrell et al., 2017; Ali et al., 2018; Munkhnasan et al., 2018; 
Rhif et al., 2020), and has been widely used in vegetation productivity estimation (Wu et al., 
2020), drought (Gong et al., 2017; He et al., 2018), desertification (Li et al., 2021), and ecological 
environment monitoring (Xie et al., 2020). Recently, scholars have studied grassland ecosystem 
changes and their driving mechanisms in different regions based on long-time series NDVI data. 
Their views can be divided into two categories: Firstly, climate change (including precipitation 
and temperature) plays the dominant role. For example, Zhao et al. (2012) found that vegetation 
degradation in the Xilin Gol grassland during 1998—2007 was significantly correlated with 
decreased precipitation. Temperature and precipitation could account for 50% of the vegetation 
variation in the Qinghai-Tibet Plateau (Wang et al., 2016). Similar results were found on the 
Loess Plateau (Naeem et al., 2021) and Ш River valley (Yan et al., 2022). Secondly, grassland 
degradation is closely related to human activities. For instance, Li et al. (2012) found that animal 
husbandry, instead of climatic factors, served as the primary vegetation degradation factor. Xilin 
Gol grassland degradation was affected by urban expansion, road construction, and mining 
activity, and that the climate change impact was less significant (Sun et al., 2017; Batunacun et al., 
2018). Especially, human activity contribution exceeded 90% of NDVI changes in the Sahel 
region's desert/steppe biome transition zone (Wu et al., 2020). Also, Pan et al. (2017) found that 
non-climatic drivers played a greater role in the grassland changes of the Qinghai-Tibet Plateau 
during 1980—2010. 

However, the abovementioned studies mostly used methods such as correlation, regression, and 
trend analysis. In fact, a strict linear relationship may not exist between vegetation growth and its 
driving factors in the complex response processes, not to mention the possible presence of 
complex interactions. For example, Shi et al. (2019) found that NDVI increased slowly overall in 
this region during 2000—2015. In areas where NDVI increased significantly, human activities 
constituted the dominant driving factor, while in areas where NDVI decreased significantly, 
precipitation changes served as the primary driving factor. To clarify the complex relationship 
between drivers and vegetation state, Wang et al. (2017) proposed a geographical detector model 
(GDM) that detects spatial heterogeneity based on spatial variance analysis, and quantified the 
relative importance of individual driving factor, and the interaction between them (Wang et al., 
2020). 

The northern Chinese grassland is ecologically sensitive, which is susceptible to climate 
change and threatened by land desertification. Natural factors and human activities dominate 
vegetation changes in this area. During the past 20 a, the average precipitation in the major sand 
and dust source regions of northern China has increased by about 20% (IPCC, 2013). Meanwhile, 
a series of ecological projects successively introduced in this area after 2000, such as Grassland 
Desertification Control, Grain for Green, and Grazing Prohibition (Hu et al., 2013). However, 
under the dual influence of climate change and human activities, the main factors driving 
vegetation change may vary and shift, which causes further uncertainty regarding the region's 
future vegetation changes. Therefore, this study focused on the Inner Mongolian Xilin Gol 
grassland, and aimed to: (1) analyze temporal and spatial variation in NDVI during 2000-2020; (2) 
identify the main driving factors of vegetation change, analyze their interactions, and determine 
the most suitable range or type of each factor in driving vegetation growth; and (3) predict the 
possible change trajectory of NDVI through multi-scenario analysis. 
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2 Materials and methods 


2.1 Study area 


Xilin Gol grassland is located in the central part of the Inner Mongolia Autonomous Region, 
China (41°35'-46°46'N, 111?09'-119*58'E) with a total area of 2.06x10? km? (Fig. 1). It has a 
mid-temperate arid and semi-arid continental monsoon climate. The area has an annual average 
temperature of 2.2?C, an average annual precipitation of about 280 mm, and an altitudinal range 
of 760-1925 m a.s.l. Its terrain is dominated by high plains, and its topography is high in the 
south and low in the north, with many low mountains and hills in the east and south, and basins 
scattered in between. The western and northern terrain is flat. Soils from the southeast to the 
northwest experience a gradual transition from chernozem to light and dark chestnut soils (Hao et 
al., 2014). The zonal vegetation is composed of grasslands, including typical meadow and desert 
grassland (Hao et al., 2014). Animal husbandry was the dominant industry, but since 2008, mining 
has emerged as the dominant economic sector, while animal husbandry has subsequently become 
the second-largest income source (Yang et al., 2011). 
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Fig. 1 Mean NDVI (normalized difference vegetation index) in the Xilin Gol grassland during 2000—2020 


2.2 Data sources 


The data used in this study mainly include NDVI, climate, topography, soil type, vegetation type, 
and human activity (Table 1). 


2.3 Data analysis 


The research framework is shown in Figure 2. Firstly, Theil-Sen (T—S) trend analysis and 
Mann-Kendall (М-К) test were performed to analyze the temporal and spatial vegetation changes. 
Secondly, anthropogenic and natural driving factors affecting NDVI changes were identified using 
GDM. Finally, multi-scenario prediction of NDVI changes was carried out using the Cellular 
Automat (CA)-Markov model. 

2.3.1 T-S trend analysis 

T-S trend analysis and M—K test can be beneficially combined to determine long-time series data 
trends, and have been successfully applied to the long-time series vegetation analysis (Meng et al., 
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Table 1 Research data used in this study 


Data type Data content Code cR Source Processing method 
. . HANTS (harmonic analysis of time series) 
Remote MODISNDVI NDVI 2000-2020 115/ и for NDVI filtering in environment for 
5 ` 8 visualizing images (ENVI) v.5.3 
The annual average temperature was queried 
Temperature (°С) ТЕМ 2000-2020 http://data.cma.cn/ by (structured query language) SQL, and Uis 
Kriging interpolation was performed in 
ArcGIS v.10.3 
The average annual precipitation was 
А DONNE T ; queried by SQL, and the Kriging 
Climate Precipitation (шт) PRE 2000-2020 http://data.cma.cn/ interpolation was performed in ArcGIS 
v.10.3 
The annual average wind velocity was 
: ; T . queried by SQL, and the Kriging 
Wind velocity (m/s) WIND 2000-2020  http://data.cma.cn/ interpolation was performed in ArcGIS 
v.10.3 
: https://ladsweb.modaps. Mosaic, masking, projection, and 
Elevation (m) BLEV 42020 eosdis.nasa.gov/search resampling in ArcGIS v.10.3 
opography spect = xtracted by aspect tool in Arc v.10. 
T hy A ASP 2020 E db l in ArcGIS v.10.3 
" Deriving slope from an elevation surface in 
Slope (°) SLP 2020 - ArcGIS v.10.3 
F : . Masking, projection, and resampling in 
Soil data Soil type SOL 2000 http://www.geodata.cn ArcGIS v.10.3 
Map of vegetation types Maski ‘acti d line i 
Vegetation Vegetation type VEG 1987 in Inner Mongolia poss ok BON) AIO а 
Autonomous Region dns 
. Masking, projection, and resampling in 
Land use type LAND 2020 http://www.resdc.cn/ ArcGIS v.10.3 
Population density POP 2000-2019 Statistical Yearbook of Joining a layer to another attribute table 
Human (person/km?) Xilin Gol League based on a common field in ArcGIS v.10.3 
activity GDP per capita GDP 2000-2019 Statistical Yearbook of Joining a layer to another attribute table 
(CNY/person) Xilin Gol League based on a common field in ArcGIS v.10.3 


Statistical Yearbook of 
Xilin Gol League 


Livestock number 


(head) LIV 2000-2019 


Note: - means no data. 
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Fig. 2 Research framework in this study 
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2020; Sun et al., 2021). T-S trend analysis, a non-parametric statistical method effectively avoids 
outliers and measurement errors without being interfered by outliers. It is therefore suitable for time 
series trend analysis with outliers (Wu et al., 2015). The calculation equation is as follows: 


X;—X 
ове = median[ i } 1<i<j<n, (1) 


j-i 
where x; and x; are the time series values in year i and year j, respectively; n is the time series data 
length (21 a in this case); median is the time series median; and Swpvi is the NDVI time series 
variation. When Snpvi>0, NDVI increases; and when Snpy1<0, NDVI decreases. The M-K test (а 
nonparametric test) is used to test the variation of series and the occurrence time of sudden 
changes (Mann, 1945; Kendall, 1990). 

Referring to the classification methods used by Yang et al. (2019) and Jiang et al. (2015), this 
study introduced three area types, i.e., improved area (Swpvr 0.0005), unchanged area 
(—0.0005<Snpvi<0.0005), and degraded area (Snpvi< —0.0005). The significance level (a=0.05) 
was tested by M-K test, and the results were divided into significant (|Z|>1.64) and 
non-significant changes (|7|<1.64). 

2.3.2 GDM 

GDM is a Statistical approach that reveals the drivers behind factors by detecting their 
hierarchical spatial heterogeneity (Chen et al., 2020; Liu et al., 2020), and was developed by 
Wang et al. (2020). This model includes factor, interaction, risk, and ecological detectors. The 
first three detectors were applied in this study, including factor detector, interaction detector, and 
risk detector. Factor detector was used to detect dependent variable spatial heterogeneity (NDVI 
in this study), and to analyze the influence of each driving factor on NDVI spatial distribution 
through a q value comparison: 


L 
Е У Мо _ 88W 
d 2 


* h =1, PA vta L, (2) 

No SST 
SSW = Y, Мной, (3) 
SST = No?, (4) 


where L is the number of classes of dependent variable Y or factor X; № and N are the numbers of 
units in class h and the entire area, respectively; о?» and o? are the variances of Y in class h and 
the entire area, respectively; SSW and SST are the intra-group sum of squares and the total sum 
of squares, respectively. The value g is the degree to which the spatial heterogeneity of vegetation 
change Y can be explained by factor X, and has a range of 0—1. A greater q value means that the 
driving factor has a greater influence on vegetation change (Zheng et al., 2021). In this study, the 
q value indicated the spatial distribution consistency between NDVI change and their drivers. 

Interaction detector was used to identify the interaction between every two factors (independent 
variables), and to evaluate the explanatory power of human activities and natural factors, jointly 
(enhanced or weakened) or separately, in explaining the spatial NDVI distribution. The 
determination method is shown in Table 2. 


Table2 Interaction relationship of each factor 


Description Interaction Description Interaction 
q(X1 0 X2)«min(q(X1), q(X2)) Weaken, nonlinear q(X1 N X2)=q(X1)+q(X2) Independent 
min(q(X1), q(X2)«q(X1 N X2)<max(q(X1), g(X2)) Weaken, univariate q(X1 0 X2)>q(X1)+q(X2) Enhance, nonlinear 
q(X1 0 X2)*max(q(X1), q(X2)) Enhance, bivariate 


Note: ХІ and X2 are the driving factors of vegetation degradation. The symbol " N "denotes the interaction between ХІ and X2. 


Risk detector was used to identify suitable and unsuitable driving factor ranges to NDVI 
(dependent variable). A test was utilized to determine whether the mean NDVI value would vary 
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across different zones for the same factor (Wang et al., 2019). 

Analysis scale differences may affect GDM results to some extent. In this study, 12 driving 
factors were selected, namely, temperature, precipitation, wind velocity, elevation, aspect, slope, 
soil type, vegetation type, land use type, population density, GDP per capita, and livestock 
number (Table 1; Fig. 3). The highest driving factor g value served as the criterion for 
establishing the optimal spatial scale (Fu et al., 2018; Table S1). After scale comparison, a value 
of 7000 m was the best spatial scale for the hierarchical spatial heterogeneity analysis (Fig. S1). 
2.3.3 Simulation and prediction of NDVI change 
(1) CA-Markov model 

CA-Markov model effectively combines the advantages of CA and Markov models. CA 
usually includes four basic elements, i.e., unit, state, adjacent range, and transformation rule 
(Wang et al., 2018). It can be expressed as follows: 


S" = fys, (5) 
where S is the finite set of discrete cell states; f is the transfer function that defines the cell 
transition from time f to time 1+1; and N is the neighborhood of the cell. 

A Markov model is a special stochastic motion process with non-aftereffect property and 
stability (Wang et al., 2018). According to Wang et al. (2018), we used a state transition matrix to 


simulate the dynamic future NDVI change (Culik et al., 1990; Wang et al., 2018; Aguejdad, 2021). 
Its expression is as follows: 
[Pi Р, Fin 
P Py ... P 
P=|P |= 21 22 2m , 0x P, <1, УРЕ (6) 
| Pat Ри Eel Pam 


where P is the state transition matrix; and P; is the transition probability from i^^ level to j® level 
in one cycle. According to the study of Wang et al. (2018), S, is the NDVI distribution of a given 
pixel at the initial time, and the cycle transfer probability matrix is obtained by multiplying with 
P".'The NDVI distribution after п cycles then becomes: 


S,=S,P". (7) 


CA-Markov model is implemented by GeoSOS-FLUS software (http://www.geosimulation.cn/ 
flus.html), which is based on the suitability probability calculation module of neural networks, 
and can quickly obtain the suitability probability of various land distribution types. 

(2) NDVI classification 

Since raster data with a discrete space and state must be used as input for CA—Markov model, 
we divided NDVI data into five levels according to the classification standard proposed by Wang 
et al. (2018): low vegetation coverage area (NDVIxO0.3), relatively low vegetation coverage area 
(0.3<NDVIS0.4), medium vegetation coverage area (0.4<NDVISO0.5), relatively high vegetation 
coverage area (0.5<NDVIS0.6), and high vegetation coverage area (NDVI>0.6), with assignments 
of 1—5, respectively. 

(3) Modeling and multi-scenario settings 

Taking NDVI spatial distribution maps in 2010 and 2015 as the basic data, this study screened 
the driving factors based on GDM analysis results, and calculated NDVI transition matrix and 
suitability map using GeoSOS-FLUS software. According to transformation rule, we simulate 
NDVI distribution in 2020 by CA—Markov model, and used Kappa precision and Figure of Merit 
(FOM) (the smaller the value, the higher the simulation accuracy) analysis to compare simulated 
values with actual values, thus testing model reliability in terms of dynamic vegetation change 
prediction. Taking 2030 as the target year, the created CA—Markov model was used to design the 
following seven development scenarios by adjusting the range size of dominant driving factors 
(Table 3). 
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Table3 Parameter settings in different scenarios 


Scenario Name Code Parameter and planning policy 
Using the seven driving factors screened in Figure 
6, keeping the model parameter settings 
Business as usual Business as usual BAU unchanged, providing no processing for the 
driving factors, and following the  inertial 
development law of NDVI. 
: : Increasing the overall wind velocity value by 
оао wae 25%, while keeping other factors constant. 
Reduced wind velocity WIN- Reducing the overall wind velocity value by 25%, 
while keeping other factors constant. 
Scenario of climate change 
RE Increasing the overall precipitation value by 25%, 
Increased precipitation PRES while keeping other factors constant. 
-— Reducing the overall precipitation value by 25%, 
Она RRE- while keeping other factors constant. 
Scenario of economic priority Increased livestock number PD+ Increasing the: value: of livestock. by:507%, While 
keeping other factors constant. 
Scenario of ecological protection Reduced livestock number PD- Reducing the value НОО while 


(a) 


Annual mean 
temperature (°C) 


E 5.82 
1.98 


Elevation (m) 


- 1893 
763 


Soil type 


=æ Lake 


(0) 


(person/hm?) 
- 44.7 
1.0 


0 200 km 
— 


(b) 


keeping other factors constant. 


(c) 


= Gray forest soil 
m Chernozem 

æ Chestnut soil 

= Brown soil 

ив Gray meadow soil 
= Boggy soil 

= Solonchak 

æ Aeolian sandy soil 


Population density 


Mean annual 
precipitation (mm) 


| И 
139.23 


Aspect 

m Flat 

æ North 

æ Northest 
æ East 

m Southest 
= South 

m Southwest 
= West 

= Northwest 
за North 


Vegetation type 

m Forest 

=æ Bushwood 

c3 Grassland 

= Desert 

= Low wetland 
vegetation 

™ Artificial vegetation 

æ Unvegetation 


(k) (1) 


GDP per capita 
(CNY/person) 


= 94,753 
24,176 


N 


A 


Mean wind 
velocity (m/s) 


= 3.92 
2.48 


Slope (°) 
26.96 


0.00 


Land use type 

= Crop land 

™ Forest land 

æ Grassland 

m Water body 

™ Urban and built-ups 
= Unused land 


M Livestock 
number (head) 
8,939,494 


lm 35,749 


Fig. 3 Map of driving factors of NDVI change in the Xilin Gol grassland. (a), annual mean temperature; (b), 
mean annual precipitation; (c), mean wind velocity; (d), elevation; (e), aspect; (f), slope; (g), soil type; (h), 
vegetation type; (i), land use type; (j), population density; (k), GDP per capita; (1), livestock number. 
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3 Results 


3.1 Spatiotemporal dynamic changes in NDVI 


Inter-annual NDVI varaitaion of the Xilin Gol grassland fluctuated between 0.25 and 0.35, 
showing an overall increasing trend (Fig. 4). The minimum value occurred in 2001, while the 
maximum values occurred in 2016. As seen from inter-annual changes, NDVI gradually increased 
after 2000, and fluctuated greatly during 2007—2013 and 2015—2018. 


0.36 
0.35 
0.34 


y=0.0027x+0.2898 


NDVI 
> 
= 


R-0.76 


0.29 
0.28 
0.27 


0.26 
2000 2002 2004 2006 2008 2010 2012 2014 2016 2018 2020 


Year 
Fig.4 Inter-annual NDVI variation of the Xilin Gol grassland during 2000—2020 


NDVI variation presented a spatially heterogeneous distribution. During 2000—2020, the areas 
with improved vegetation coverage accounted for 84.5396, while the degraded areas accounted for 
7.43% (Table 4). The areas with significantly improved vegetation were distributed in the 
southeastern and northeastern regions of the Xilin Gol grassland, as well as some areas of central 
region. The areas with severely degraded vegetation were scattered in the eastern, northern, and 
central regions of the Xilin Gol grassland (Fig. 5). 


Table4 Statistics of NDVI trend 


Swpvi Z value Trend of NDVI Area percentage (%) 
S< —0.0005 Zs —1.64 Significant decrease 0.67 
S< —0.0005 —1.64<7<1.64 Slight decrease 6.76 
—0.0005—0.0005 —1.64<7<1.64 Stability 8.04 
$>0.0005 —1.64<7<1.64 Slight increase 58.11 
5>0.0005 721.64 Significant increase 26.42 
z. 115°Е 120°Е 


45N  46?N 
45?N  46?N 


44?N 


44°N 


Trend of NDVI 

m= Significant decreasing 
= Slight decreasing 

© Stability 

= Slight increasing 

= Significant increasing 


43?N 
42°м 43?N 


42?N 


115?B 120°E 
Fig.5 Trend of NDVI of the Xilin Gol grassland during 2000—2020 
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3.2 Drivers of vegetation change 


3.2. Effect of single driving factor on vegetation change 

Factor detector revealed the driving effects of natural and human activity factors on NDVI change. 
These single factors were ranked as follows in a descending order: wind velocity? precipitation? 
soil type>livestock number^temperature» population density» GDP per capita? vegetation type? 
slope>land use type» DEM? aspect (Fig. 6). Among them, the effects of wind velocity (920.6236) 
and precipitation (420.6190) оп NDVI exceeded 0.6, and were therefore the dominant factors 
affecting vegetation status of the Xilin Gol grassland. The explanatory power of livestock number 
in explaining NDVI change in the Xilin Gol grassland exceeded 0.4 (920.4490), suggesting that 
grazing is also an important factor for NDVI change. 


Driving factor 


0.3327 


0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 
q value 


EB Human activity driver factor Natural driver factor 
Fig. 6 Single factor affecting vegetation change in the Xilin Gol grassland during 2000—2020 and its q value. 
TEM, temperature; PRE, precipitation; WIND, wind velocity; ELEV, elevation; ASP, aspect; SLP, slope; SOL, 
soil type; VEG, vegetation type; LAND, land use type; POP, population density; GDP, GDP per captital; LIV, 
livestock number. The abbreviations are the same in Figure 7. 


3.2.2 Effects of interactions among driving factors on vegetation change 


Interaction detector evaluated the interaction effects among driving factors on NDVI change 
(Table S2; Fig. 7). First, the interaction q values among driving factors all exceeded that of single 
factor (Table S2), showing mutual or nonlinear enhancement. Among natural factors, average 
annual precipitation and wind velocity had the strongest interaction (g=0.7456). Secondly, the 
interaction explanatory power between average annual precipitation and temperature in 
explaining the spatial NDVI distribution (4=0.7192) also exceeded 0.7. Significant interactions 
also existed among factors such as wind velocity, temperature, precipitation, and soil use type 
(a=0.05, with g values above 0.5000 in each case). The interaction g values between human 
activity and other factors all showed an increasing trend, and the strongest interaction was 
observed between average annual precipitation and livestock number (g=0.7182). In this study, 
after two-factor interactions, nonlinear enhancement was uniformly found in the interactions 
between slope and other factors (Fig. 7). This result indicated that slope, as an indirect driver of 
vegetation growth, affects NDVI by affecting other factors. 

3.2.3 Suitability analysis on the influencing factors of NDVI 

Risk detector can be used to detect the suitable range or type of each driving factor in driving 
NDVI increase or decrease. The greater the mean NDVI value of different zones, the more 
favorable the range for vegetation growth. Detection results showed that precipitation (suitable 
range of 303.89—393.36 mm), elevation (suitable range of 1256—1695 т), slope (suitable range of 
8.22°—16.18°), and population density (suitable range of 3.6—13.9 person/kn?) were all positively 
correlated with NDVI (Fig. 8; Table 5), while temperature (suitable range of 2.74°C—3.33°C) and 
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Fig. 7 Interactions among driving factors of NDVI 
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Fig. 8 Risk detection results of driving factors of mean NDVI. (a), temperature; (b), precipitation; (c), wind 
velocity; (d), elevation; (e), aspect; (f), slope; (g), soil type; (h), vegetation type; (i), land use type; (j), population 


density; (k), GDP per capita; (1), livestock number. 


Table5 Suitable range/type of each driving factor (9596 confidence level) 


Factor Indicator Suitable range/type Mean NDVI 
TEM Temperature (°C) 2.74—3.33 0.2262 
PRE Precipitation (mm) 303.89—393.36 0.2540 

WIND Wind velocity (m/s) 2.48—2.84 0.2487 

ЕГЕУ Elevation (т) 1256-1695 0.2284 
ASP Aspect North 0.1994 
SLP Slope (°) 8.22-16.18 0.3093 
SOL Soil type Gray forest soil 0.3253 
VEG Vegetation type Broad-leaved forest 0.3321 

LAND Land use type Forest land 0.3127 
POP Population density (person/km?) 3.6-13.9 0.2262 
GDP GDP per capita (CNY/person) 47,431—74,442 0.2288 
LIV Livestock number (head) 1,486,331-2,570,527 0.2309 
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wind velocity (suitable range of 2.48—2.84 m/s) exhibited significant negative correlations with 
NDVI. The vegetation, soil, and land use types most suitable for vegetation growth were 
deciduous broad-leaved forest, gray forest soil, and forest land, respectively. 


3.3 Simulation of spatiotemporal change of NDVI 


3.3.4 Simulation process and precision 

Based on mean NDVI distribution map in 2010 and 2015 (Fig. 9a and b), we used the 
GeoSOS-FLUS v.2.3 software to generate NDVI transition matrix and transition probability 
(Tables S3 and S4). According to Figure 6, we selected driving factors with an explanatory power 
exceeding 0.3 for NDVI changes (temperature, precipitation, wind velocity, soil type, population 
density, GDP per capita, and livestock number) for predicting NDVI changes. NDVI 
transformation suitability probability map was created, and NDVI spatial distribution map in 
2020 was acquired through simulation (Fig. 9d). A comparison between actual and simulated 
mean NDVI spatial distribution maps in 2020 (Fig. 9c and d) resulted in a Kappa coefficient of 
0.7023 and an FoM coefficient of 0.1526, which suggested a high simulation precision. 


Pz 


Real spatial 
distribution of 
NDVI in 2010 


Real spatial 
distribution of 
NDVI in 2015 


Simulated spatial distribution 

of NDVI in 2020 

ви Low vegetation coverage 
area (NDVI<0.3) 
Relatively low vegetation 
coverage area (0.3<NDVI<0.4) 
Medium vegetation coverage 


© 

Real spatial E 
area (0.4<NDVI<0.5) 

EH 

m 


distribution of 
NDVI in 2020 Relatively high vegetation 
coverage area (0.5<NDVI<0.6) 
0 200 km à High vegetation coverage 

area (NDVI>0.6) 


Fig.9 Real spatial distribution of NDVI in 2010 (a), 2015 (b), 2020 (c) and its simulated result in 2020 (d) 


3.3.2 Simulated spatial distribution of NDVI under different scenarios 


Figure 10 showed simulated spatial change of NDVI under different scenarios in 2030. NDVI 
transition path was low vegetation coverage area-relatively low vegetation coverage area— 
medium vegetation coverage area-relatively high vegetation coverage area-high vegetation 
coverage area under the scenarios of reduced wind velocity (WIN-, Fig. 10c), increased 
precipitation (PRE+, Fig. 10d), and ecological protection (PD-, Fig. 10g). Compared with the 
business as usual (BAU, Fig. 10a), the above scenarios showed the increases in high vegetation 
coverage area by 18,195, 11,206, and 16,174 km?, respectively (Fig. 11). Overall, NDVI showed 
an increasing trend (vegetation restoration), and the increased areas located in the northeast (West 
Ujimqin Banner and East Ujimqin Banner) and southeast (Duolun County and Taipusi Banner) of 
the Xilin Gol grassland. 

Under the scenarios of increased wind velocity (WIN+, Fig. 10b), decreased precipitation 
(PRE-, Fig. 10e), and economic priority (PD+, Fig. 10f), NDVI transition path showed an 
opposite trend. Compared with the business as usual (BAU, Fig. 10a), these scenarios showed 
decreases in high vegetation coverage area by 1000, 854, and 783 km?, respectively, and 
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decreases in relatively high vegetation coverage area by 1895, 2630, and 1938 km/, respectively 
(Fig. 11). Overall, NDVI showed a decreasing trend, suggesting that vegetation coverage was 
significantly reduced under climate conditions that were not conducive to vegetation growth. 
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Fig. 10 Simulated spatial distribution of NDVI in the Xilin Gol grassland under different scenarios in 2030. (a), 
BAU; (b), WIN+; (c), WIN-; (d), PRE+; (e), PRE-; (f), PD+; (g), PD-. The detailed scenarios are presented in 
Table 3. 
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Fig. 11 Area change of NDVI at each level under different scenarios (WIN+, WIN-, PRE+, PRE-, PD+, and PD-) 
in 2030. The detailed scenarios are presented in Table 3. 


4 Discussion 


4.4 Dynamic changes in NDVI 


T-S trend analysis found that, since 2000, vegetation in the Xilin Gol grassland has been in a 
state of overall recovery and partial degradation, and that vegetation changes have shown obvious 
spatial heterogeneity, which is consistent with the results of Batunacun et al. (2018) and Shi et al. 
(2019). Variation in NDVI demonstrated a distribution pattern of high in the northeast and low in 
the southwest, and the areas where the vegetation coverage has increased significantly locate 
mainly in the northeastern (East Ujumqin Banner), central (Xilinhot City and Abag Banner), and 
southeastern regions of the Xilin Gol grassland (Duolun County, Taipusi Banner, and Zhenglan 
Banner). Although the Xilin Gol grassland has been restored to a certain extent during the past 21 
a, serious degradation has still occurred in some areas. These areas mainly located in West 
Ujimqin Banner and East Ujimqin Banner, and scattered in the Xilinhot City and Sonid Left 
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Banner. Grassland degradation in these areas is associated with local land use changes, especially 
the development of mining industry (Hu et al., 2013). For example, during open-pit mine (coal 
mines and quarries) exploitation, wastage is stripped off in large quantities during the mining 
process to form numerous dumps, which occupies grasslands and destroys vegetation, resulting in 
a degraded ecological environment (Dai et al., 2014). Moreover, rapid development of industrial 
infrastructure related to mining has further aggravated grassland fragmentation, eventually 
leading to grassland degradation (Qian et al., 2014). 


4.2 Driving mechanisms of NDVI changes 


Factor detection revealed that precipitation and wind velocity were the main natural factors 
affecting vegetation changes, while temperature did not play a significant role. This finding is in 
accordance with Zhang et al. (2006), Zhao et al. (2012), Wang et al. (2016), and Naeem et al. 
(2021). Precipitation is positively correlated with NDVI (Zhao et al., 2012). Precipitation is the 
leading factor of vegetation growth and productivity improvement (Liu et al., 2019), and 
determines the distribution space and growth state of vegetation. The warming and humidification 
change of climate in northern China in recent decades (IPCC, 2013) has provided a more 
favorable growth environment for vegetation in the Xilin Gol grassland. Among them, Duolun 
County, Taibus Banner and the northeast of East Ujimqin Banner have the most significant 
improvement. A wind velocity increase leads to increased soil evapotranspiration. Consequently, 
decreased soil water content induces water stress, which is detrimental to net primary grassland 
productivity, and ultimately affects grassland NDVI (Gardiner et al., 2016). Excessive wind 
velocity blows soils away and damages plant seedlings, resulting in decreased land productivity 
and serious ecological problems (Wang et al., 2016; Zou et al., 2018). In contrast, topographic 
factors and vegetation types affect vegetation changes to lesser extents (with q values less than 
0.2000). 

Human activities also affected the spatial distribution of NDVI. A very strong positive 
correlation existed between grassland degradation and livestock number (Xie et al., 2012; Hao et 
al. 2014), and residual trend analysis showed that animal husbandry was the main driver of 
vegetation change in the Xilin Gol grassland during 1981—2006 (Li et al., 2012). In this regard, 
Sun et al. (2017) and Batunacun et al. (2018) reported similar results. This study also concluded 
this, that is, livestock number (920.4490) and population density (420.3220) are the main human 
factors affecting vegetation change. Nonetheless, risk detection indicated that NDVI was 
positively correlated with population density, and that local large-scale vegetation restoration 
(increased NDVI) might be related to ecological restoration projects and policies, intending to 
balance grasslands and livestock (Dong et al., 2007; Shi et al., 2019). Particularly, a series of 
ecological restoration projects, such as Grassland Desertification Control, Grain for Green, 
Ecological Immigration Policy, and Beijing-Tianjin Sandstorm Source Control (Hu et al., 2013) 
had promoted local vegetation restoration (Shi et al., 2019) since 2000. Additionally, grazing has 
been prohibited in severely degraded grasslands by executing a policy of grazing rest for 
moderately degraded grasslands, and by implementing rotational grazing for mildly degraded 
grasslands since 2002 (Dong et al., 2007; Shi et al., 2019). The implementation of the 
above-mentioned projects and policies effectively reduced the pressure imposed by human 
activities on NDVI changes. Consequently, the grassland degradation drivers gradually shifted 
from anthropogenic (such as grazing) to natural (climate change) factors. Correspondingly, the 
effects of human activities were weakened, while those of climatic factors were strengthened. 

NDVI changes in this area were caused by the superposition of multiple factors. The two-factor 
superposition of driving factors increased their explanatory power in the spatial distribution of 
NDVI (Fig. 7); that is, the interactions between climatic factors, as well as between climatic 
factors and other factors, enhanced their explanatory power. For example, the g values of average 
annual wind velocity and precipitation were both «0.6500 in explaining NDVI separately (Fig. 6). 
However, the strongest interaction occurred between them (4=0.7456), followed in succession by 
the interaction between average annual precipitation and annual average temperature (4=0.7192). 
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The combined action of multiple factors ultimately affects vegetation growth (Hao et al., 2017; 
Qiao et al., 2018). In this study, soil type combined with other factors also exerted important 
effects on vegetation growth and rainwater re-use (Piao et al., 2006; Otgonbayar et al., 2017). 


4.3 Policy implications 


The factors driving the changes in grassland ecosystem services are distinctively different at 
different time scales. Policy makers should therefore formulate suitable and time-specific 
measures for grassland ecosystem management (Hao et al., 2017). Although the Xilin Gol 
grassland is progressing towards recovery, its vegetation status is still threatened by climate 
change and human activities, mainly due to its fragile ecological environment. The multi-scenario 
simulation results provide reference and data support for designing grassland restoration policies 
to coordinate regional socio-economic development and ecological governance. For example, 
under the scenarios of reduced wind velocity, increased precipitation, and reduced livestock 
population, overall vegetation recovery exists. In contrast, under the scenarios of increased wind 
velocity, reduced precipitation, and increased livestock population, significant vegetation 
degradation occurs. 

Grasslands are a renewable resource. Under reasonable management, grasslands can be 
continuously regenerated and reproduced, and be put into sustainable use by human beings. If this 
is not done correctly, such grasslands will degrade or even disappear. Thus, for future ecological 
management, governments should prioritize the determination of grassland areas with high water 
availability, and correlate grass availability with livestock population density. Together with this, 
grassland management and animal husbandry should be deployed scientifically. Depending on the 
specific water resource availability in the study area, grassland carrying capacity should be 
determined, and the avoidance of livestock overload should be taken as the primary measure to 
curb grassland degradation. Simultaneously, animal husbandry development should be adopted as 
the breakthrough point for agriculture and animal husbandry industrialization, and artificial grass 
planting should be vigorously implemented to promote such activities and to balance grassland 
and livestock biomass. 

Detailed recommendations are as follows: according to the variation range of annual 
precipitation and wind speed, the range of grassland carrying capacity should be comprehensively 
evaluated. At the same time, according to the natural and socio-economic conditions of different 
regions, the restoration and management of severely degraded grassland should be focused on, 
and the ecological protection of moderately and mildly degraded grassland should be 
appropriately strengthened to realize the sustainable development of agriculture and animal 
husbandry. 


4.4 Limitations of the study 


This study utilized the MODIS NDVI dataset to analyze the distribution of, and changes in 
vegetation of the Xilin Gol grassland during 2000—2020. The spatial resolution of this dataset is 
relatively low, which means that it may not be an ideal choice for some areas, and may end up 
exaggerating the research results of these areas. For example, in areas with dense vegetation 
coverage, the data quality was acceptable, while in areas with sparse vegetation coverage, NDVI 
was greatly affected by soil. Therefore, NDVI data of Landsat or Sentinel-2 can be considered in 
the future. Additionally, due to the rapid fluctuation of annual plants with precipitation, using 
NDVI as an indicator of vegetation status may produce erroneous estimates of vegetation 
restoration status in areas with low vegetation coverage (Tao et al., 2020). Therefore, net primary 
productivity and enhanced vegetation index can be introduced to improve future research. 

The selection of driving factors did not consider the influence of groundwater, soil water 
content, market, or policy factors. In the future, more remote sensing climatic data (such as 
groundwater dynamic changes, soil moisture, evapotranspiration, and drought index) and policy 
factors can be taken into account for an in-depth analysis of NDVI spatial distribution drivers. 

The CA-Markov model could predict the changes in vegetation in the Xilin Gol grassland, but 
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due to the interactions among data, models, driving factors, and other aspects, there were errors in 
prediction results. Limited by the model transformation rule, it was necessary to convert 
continuous NDVI data into five levels (discrete type), which led to precision loss in the data input 
into the model, affecting simulation results. Subsequent research can incorporate long short-term 
memory (LSTM) into the Markov transformation rule, and repeatedly consider the phase of 
driving factors to improve model precision. 


5 Conclusions 


This study revealed the temporal and spatial distribution characteristics of NDVI in the Xilin Gol 
grassland during 2000—2020, its natural and human activity drivers, and the effects of interactions 
between drivers on NDVI changes. The effects of dominant driving factors on the spatial 
distribution of NDVI under multiple scenarios were also simulated using CA—Markov model. 
Vegetation of the Xilin Gol grassland had an overall improvement, however, it was still being 
degraded locally. Temporally, NDVI in the Xilin Gol grassland generally showed a gradual 
increase during 2000—2020. Spatially, NDVI exhibited a distribution pattern of high in the 
northeast and low in the southwest. After 2000, natural factors were the dominant factors 
affecting vegetation coverage in the Xilin Gol grassland, but the influence of human activities 
could not be ignored. Coordinating the relationship between human activities and natural 
environment is greatly significant to vegetation restoration and ecological protection. 
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Appendix 
Table S1 Zoning effect of geographical detector 
Category Factor Method qvalue | Category Factor Method q value 
Natural break 0.3333 Natural break 0.0981 
Temperature Quantile 0.3301 Elevation Quantile 0.111 
Geometrical Geometrical 
| 0.2884 . 0.0781 
interval interval 
Topography 
Natural break 0.5995 Natural break 0.1545 
7 Quantile 0.6262 Quantile 0.1213 
Precipitation Slope 
Geometrical Geometrical 
ў 0.5770 . 0.1232 
interval interval 
Climate Natural break 0.6351 Natural break 0.2305 
. Quantile 0.6447 Population Quantile 0.3188 
Drought index : 
Geometrical density Geometrical 
. 0.5899 | 0.3177 
interval interval 
Natural break 0.6283 Natural break 0.1256 
. Quantile 0.6313 JAnthropogenic GDP Quantile 0.3242 
Wind speed cms : 
Geometrical owe per capita Geometrical 
б 0.6162 . 0.3060 
interval interval 
Natural break 0.1583 Natural break 0.4491 
Soil pH value Quantile 0.1389 Livestock Quantile 0.2921 
Geometrical Geometrical 
: 0.1595 ў 0.3121 
interval interval 
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Fig. 51 Scale effect of geographic detector result (q value and 90% quantile of q value). х1—х14 are the detailed 


driving factors. 
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Table S2 Interaction of driver factor 
C=q(X1NX2) A=q(X1) B=q(X2) Conclusion Interpretation 
х10х2=0.7192 0.3327 0.6190 С<А+В:С>А.В * 
х10х3=0.7094 0.3327 0.6379 С<А+В;С>А,В * 
х10х4=0.6955 0.3327 0.6236 С<А+В;С>А,В E 
х10х5=0.4284 0.3327 0.0987 С<А+В;С>А,В z 
x1Nx6=0.3469 0.3327 0.0051 C>A+B;C>A,B t 
х10х7=0.4776 0.3327 0.1549 С<А+В;С>А,В i 
х10х8=0.5781 0.3327 0.4635 С<А+В;С>А,В * 
х1Пх9=0.4402 0.3327 0.1653 C<A+B;C>A,B * 
х10х10=0.5063 0.3327 0.1830 С<А+В;С>А,В * 
х10х11=0.4832 0.3327 0.1428 С<А+В;С>А,В 
х10х12=0.6069 0.3327 0.3220 С<А+В;С>А,В i 
x1Nx13=0.6562 0.3327 0.3217 C<A+B;C>A,B е 
х10х14=0.6323 0.3327 0.4490 С<А+В;С>А,В * 
х2Пх3=0.6586 0.6190 0.6379 С<А+В;С>А,В * 
х2Пх4=0.7456 0.6190 0.6236 C<A+B;C>A,B 
x2x5=0.665 1 0.6190 0.0987 C<A+B;C>A,B = 
х201х6=0.6285 0.6190 0.0051 С>А+В;С>А,В T 
х201х7=0.6529 0.6190 0.1549 С<А+В;С>А,В m 
х2Пх8=0.7063 0.6190 0.4635 C<A+B;C>A,B * 
х2Пх9=0.6863 0.6190 0.1653 C<A+B;C>A,B * 
х2Пх10=0.6597 0.6190 0.1830 С<А+В;С>А,В + 
х2Пх11=0.6699 0.6190 0.1428 C<A+B;C>A,B и 
х20х12=0.6754 0.6190 0.3220 С<А+В;С>А,В * 
х20х13=0.7133 0.6190 0.3217 С<А+В;С>А,В * 
х2Пх14=0.7182 0.6190 0.4490 C<A+B;C>A,B * 
x3x4=0.7179 0.6379 0.6236 C<A+B;C>A,B 
х3Пх5=0.6690 0.6379 0.0987 C<A+B;C>A,B * 
x3x6=0.6488 0.6379 0.0051 C>A+B;C>A,B d 
x3x7=0.6856 0.6379 0.1549 C<A+B;C>A,B б 
x3Nx8=0.7272 0.6379 0.4635 C<A+B;C>A,B * 
х3Пх9=0.6919 0.6379 0.1653 С<А+В;С>А,В * 
х3Пх10=0.6936 0.6379 0.1830 C<A+B;C>A,B * 
x3Mx11=0.7019 0.6379 0.1428 C<A+B;C>A,B i 
x3Mx12=0.6754 0.6379 0.3220 С<А+В;С>А,В * 
х3Пх13=0.6878 0.6379 0.3217 C<A+B;C>A,B е 
х3 х14=0.6909 0.6379 0.4490 С<А+В;С>А,В * 
х40х5=0.7112 0.6236 0.0987 С<А+В;С>А,В * 
x4(1x6—0.6343 0.6236 0.0051 С>А+В;С>А,В T 
х4Пх7=0.6696 0.6236 0.1549 C<A+B;C>A,B 2 
x4(1x8-0.7144 0.6236 0.4635 C<A+B;C>A,B " 
x4(1x9—0.6553 0.6236 0.1653 C<A+B:C>A,B Ы 
х40х10=0.6775 0.6236 0.1830 С<А+В;С>А,В * 
x4(1x11-0.6844 0.6236 0.1428 C<A+B;C>A,B * 
х40х12=0.6893 0.6236 0.3220 С<А+В;С>А,В 
х4Пх13=0.6800 0.6236 0.3217 C<A+B;C>A,B ы 
х4Пх14=0.6414 0.6236 0.4490 C<A+B;C>A,B * 
х5Пх6=0.1131 0.0987 0.0051 C>A+B;C>A,B t 
х5Пх7=0.2682 0.0987 0.1549 С<А+В;С>А,В * 
х5Пх8=0.5320 0.0987 0.4635 C<A+B;C>A,B * 
х5Пх9=0.2834 0.0987 0.1653 C<A+B;C>A,B id 
х50х10=0.2655 0.0987 0.1830 С<А+В;С>А,В Ы 
х50х11=0.2289 0.0987 0.1428 С<А+В:С>А.В * 


To be continued 


JOURNAL OF ARID LAND 


Continued. 
C=q(X1NX2) A=q(X1) B=q(X2) Conclusion Interpretation 
x5Mx12=0.3949 0.0987 0.3220 C<A+B:C>A.B * 
x5Mx13=0.5603 0.0987 0.3217 C<A+B;C>A,B * 
х5Пх14=0.5644 0.0987 0.4490 C<A+B;C>A,B Ы 
x6Mx7=0.1691 0.0051 0.1549 C>A+B;C>A,B T 
x6Mx8=0.4859 0.0051 0.4635 C>A+B;C>A,B T 
x6Mx9=0.1916 0.0051 0.1653 C>A+B;C>A,B Т 
x6Mx10=0.2029 0.0051 0.1830 C>A+B;C>A,B T 
x6Mx11=0.1590 0.0051 0.1428 С>А+В;С>А,В t 
x6Mx12=0.3342 0.0051 0.3220 C>A+B;C>A,B Т 
х6Пх13=0.3348 0.0051 0.3217 C>A+B;C>A,B T 
х6Пх14=0.4654 0.0051 0.4490 C>A+B;C>A,B Т 
x7Mx8=0.5047 0.1549 0.4635 C<A+B;C>A,B * 
xX7Mx9=0.2882 0.1549 0.1653 C<A+B;C>A,B * 
x7Mx10=0.2823 0.1549 0.1830 C<A+B;C>A,B ii 
x7Mx11=0.2424 0.1549 0.1428 С<А+В;С>А,В Ы 
x7Mx12=0.4222 0.1549 0.3220 C<A+B;C>A,B * 
x7Mx13=0.4106 0.1549 0.3217 C<A+B;C>A,B * 
x7Mx14=0.5165 0.1549 0.4490 С<А+В;С>А,В * 
x8Mx9=0.5402 0.4635 0.1653 C<A+B;C>A,B * 
x8Mx10=0.5350 0.4635 0.1830 C<A+B;C>A,B * 
x8Mx11=0.5375 0.4635 0.1428 C<A+B;C>A,B Ы 
x8Mx12=0.6301 0.4635 0.3220 C<A+B;C>A,B Ы 
x8Mx13=0.5966 0.4635 0.3217 C<A+B;C>A,B * 
х8Пх14=0.6628 0.4635 0.4490 C<A+B;C>A,B * 
x9Mx10=0.3169 0.1653 0.1830 С<А+В;С>А,В ° 
x9Mx11=0.2952 0.1653 0.1428 C<A+B;C>A,B * 
x9f1x12-0.4319 0.1653 0.3220 C<A+B;C>A,B * 
x9Mx13=0.4137 0.1653 0.3217 C<A+B;C>A,B * 
х9Пх14=0.5422 0.1653 0.4490 C<A+B;C>A,B * 
х10Пх11=0.2631 0.1830 0.1428 C<A+B;C>A,B * 
х10Пх12=0.4232 0.1830 0.3220 C<A+B;C>A,B * 
х10Пх13=0.4587 0.1830 0.3217 С<А+В;С>А,В Ы 
х10Пх14=0.5300 0.1830 0.4490 C<A+B;C>A,B Ы 
х11Пх12=0.4188 0.1428 0.3220 C<A+B;C>A,B * 
х11Пх13=0.4380 0.1428 0.3217 С<А+В;С>А,В * 
х11Пх14=0.5287 0.1428 0.4490 C<A+B;C>A,B * 
х12Пх13=0.6098 0.3220 0.3217 C<A+B;C>A,B 
х12Пх14=0.6103 0.3220 0.4490 C<A+B;C>A,B * 
х13Пх14=0.5671 0.3217 0.4490 C<A+B;C>A,B ы 
Note: * indicates that interaction of two factors is mutually reinforcing; | indicates that interaction of two factors is nonlinearly 
enhanced. C is the q value of interaction between two drivers. x1—x14 are the detailed driving factors. 


Table S3 Conversion matrix of NDVI 


NDVI type Very low Low Medium High Very high 
Very low 48,087 25,769 4551 310 31 
Low 1412 12,212 13,617 3705 326 
Medium 113 2987 15,283 12,800 2621 
High 22 342 5196 14,269 8505 


Very high 21 33 576 5160 23,856 
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Table S4 Conversion probability of NDVI 


NDVI type Very low Low Medium High Very high 
Very low 0.610647 0.327237 0.057788 0.003937 0.000391 
Low 0.045158 0.390517 0.435443 0.118467 0.010415 
Medium 0.003339 0.088354 0.452113 0.378662 0.077532 
High 0.000768 0.012079 0.183370 0.503620 0.300163 


Very high 0.000706 0.001105 0.019444 0.174058 0.804688 


